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The generation of energetic electron bunches by the interaction of a short, ultra-intense 
(/ > 10 19 W/cm 2 ) laser pulse with “grating” targets has been investigated in a regime of ultra- 
high pulse-to-prepulse contrast (10 12 ). For incidence angles close to the resonant condition for 
Surface Plasmon (SP) excitation, a strong electron emission was observed within a narrow cone 
along the target surface, with energies exceeding 10 MeV. Both the energy and the number of emit¬ 
ted electrons were strongly enhanced with respect to simple flat targets. The experimental data 
are closely reproduced by three-dimensional particle-in-cell simulations, which provide evidence for 
the generation of relativistic SPs and for their role in driving the acceleration process. Besides the 
possible applications of the scheme as a compact, ultra-short source of MeV electrons, these results 
are a step forward the development of high field plasmonics. 

PACS numbers: 


Surface plasmons (SP) are electromagnetic modes 
which can be excited at a sharp material interface hav¬ 
ing a periodic modulation, e.g. a grating, to allow phase 
matching with incident laser light. SP excitation allows 
enhancement and fine control of electromagnetic field 
coupling with structured metal surfaces, which is the ba¬ 
sis for the vibrant research field of plasmonics [l| with 
several applications ranging from light concentration be¬ 
yond diffraction limit [J to biosensors [3] and plasmonic 
chips Q. 


Extending the study of plasmonics to very high fields 
might allow novel applications, such as tapered wave¬ 
guides for high energy concentration and field amplifi¬ 
cation to extreme values jB, 6] or enhanced energy ab¬ 
sorption which could be exploited, e.g., for proton ac¬ 
celeration by integrating the gratings in complex target 
assemblies (see e.g. [3, 0)- In particular, the phase ve¬ 
locity and the longitudinal field component of SPs make 
them suitable for accelerating electrons to high energies, 
provided that the SP field is driven high enough for the 
electron dynamics to become relativistic. While a the¬ 
ory of relativistic SP is still lacking, numerical simula¬ 
tions have provided evidence of SP-related effects in the 
relativistic regime im including electron acceleration 
in grating targets at weakly relativistic intensities |ll|. 
However, most experiments so far have been restricted 
to intensities below 10 16 W/cm 2 12|-14|, unable to drive 
a relativistic SP, because of the intense prepulses inher¬ 
ent to high-power laser systems which can lead to an 
early disruption of the target structuring. The devel¬ 
opment of devices for ultra-high contrast pulses [15|, Il6 ] 
now allows to explore the interaction with solid targets 


structured on a sub-micrometric scale in the relativistic 
regime EzLIIl- In particular, an enhancement of the cut¬ 
off energy of protons accelerated from the rear surface of 
grating targets was observed and related to SP excitation 

In this Letter, we demonstrate electron acceleration 
with relativistic SPs excited on the surface of grating 
targets at laser intensities I = 5 x 10 19 W/cm 2 , cor¬ 
responding to a relativistic parameter ao — 5 where 
ao = (/A 2 /10 18 W cm _2 /im 2 ) 1 / 2 . With respect to flat 
targets, we observed a strong enhancement of both the 
energy and number of electrons emitted from gratings ir¬ 
radiated at an incidence angle close to the resonant value 
for SP excitation. Electron emission was concentrated in 
a narrow cone close to the target surface, with energies 
exceeding 10 MeV. 

The basics of SP generation and electron accelera¬ 
tion may be described as follows. At high laser in¬ 
tensities (/ > 10 18 W/cm 2 ) a solid target is ionized 
within one laser cycle, thus the interaction occurs with 
a dense plasma. Assuming a plasma dielectric function 
e(uj) = 1 — cjp/cu 2 = 1 — a (where uo p is the plasma 
frequency) the phase velocity of a SP is Vf = uj/k = 
c(a — 2) 1 / 2 /(<a — l) 1 / 2 where k is the SP wave vector 
and a > 2 holds. The condition for resonant excita¬ 
tion of a SP on a periodically modulated target (grat¬ 
ing) by an incident EM wave of the same frequency is 
A/A g = (1 — <a) 1 / 2 /(2 — a) 1 ! 2 — sin (fa) where A^ is the 
grating period and fa is the angle of incidence. Notice 
that these equations neglect thermal and collisional ef¬ 
fects as well as possible relativistic non-linearities, thus 
in principle the resonance could be expected at somewhat 
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different angles. 

Electron acceleration up to relativistic energy by the 
longitudinal field of a SP requires the phase velocity Vf to 
approach c in order to minimize dephasing, thus a 1 is 
required as expected for a solid-density plasma and opti¬ 
cal frequencies. The basic process may be described very 
similarly to the well-known acceleration in wake plasma 
waves [ 2 QJ, but with the difference that the transverse 
field component of the SP pushes electron on the vacuum 
side, so that the process is two-dimensional and eventu¬ 
ally the electrons are emitted at some angle with respect 
to the SP propagation direction. In a frame L' moving 
with velocity vj = Vf y with respect to the laboratory 
frame L, the SP field is electrostatic in the vacuum re¬ 
gion (x < 0 ) and can be derived from the potential 


$' = -^ e k ' x sin k'y' , (1) 

k' 

where 7 / = (1 — v 2 /c 2 ) -1 / 2 = (a — l) 1 / 2 , k' = kj 7 , 
and E s P is the amplitude of the longitudinal SP field 
(E y ) in L. Thus, in the L' frame the process is simply 
described as the electron going downhill the potential 
energy —e3>'. Due to the evanescent field component 
E' x = — d x <$>' , electrons are predominantly accelerated 
towards the vacuum side with velocity almost normal to 
the x = 0 surface. The condition of “optimal” injection 
corresponds to an electron placed initially in L' at the 
top of the potential hill (V = 0 , y' = n/2k') with v' y = 
0, i.e. with an initial velocity vj in L. Such electron 
will acquire in L' the energy W' = eE SP ^f/k. If a SP = 
eE SP /m e ujc > 1 then W' m e c 2 . In this limit, the 
energy-momentum in L' is p f ~ (IT', FF'/c, 0, 0) and thus 
pn ~ ( 7 /W 7 , IT'/c, IfW'/c, 0) in L. The final energy and 
emission angle a e are given by 


eEs-p'j 2 

of ~ ~ m e c a SP a , 


tan0 e = — ~ 7 f 1 . (2) 
Py 1 


Thus, strongly relativistic electrons ( £f m e c 2 ) are 

emitted at small angles 0 e , i.e. close to the target surface. 
The acceleration length £ a = £f /eE SP ~ \a/2ir , showing 
that electrons may reach the highest energy over a few 
microns. 

The experiment was carried out at the CEA Saclay 
Laser Interaction Center (SLIC) facility with the UHI100 
laser system (see [ 21 ] for a preliminary presentation of 
the experimental results). UHI100 provides a 25 fs laser 
pulse with a peak power of 100 TW and a wavelength 
A ^ 0.8 pm. A double plasma mirror [Hj yielded a very 
high pulse contrast (^ 10 12 ). The pulse was focused on 
target with an off-axis f/3.75 parabola in a focal spot 
of ^ 4 pm FWHM containing ~ 60% of the total en¬ 
ergy in the 1 /e 2 spot diameter, which lead to an average 
intensity of ^ 5 x 10 19 W/cm 2 . Focal spot optimisa¬ 
tion was performed with an adaptive optical system. P- 
polarization was used throughout the experiment. The 



FIG. 1: Schematic view of the experimental setup. The 2D 
top view (left) shows the position of the diagnostics. The 3D 
sketch (right) shows the adopted convention for the angles (j) 
and 0. 


schematic view of the experimental setup is shown in 
Figffl A compact CMOS based spectrometer, specifi¬ 
cally designed for this experiment, was mounted on a 
motorised tray able to change the angle </> spec within the 
range 0 ° — 60° while remaining aligned to the interac¬ 
tion center. The entrance lead pinhole had a diameter 
of 500 /im and was placed at 8 cm from the interaction 
point. A pair of permanent magnets dispersed the elec¬ 
trons on the large (49.2 x 76.8 mm 2 ) triggered 12bit 
CMOS with 48 pm pixel size. The energy detection 
range was ^2 — 30 MeV. A scintillating Lanex screen 
(16 x 5 cm 2 ) was used to collect the spatial distribution 
of the electrons in the full angular range </> = 0° — 90°. 
The screen was placed with an angle of 45° with respect 
to the target and its center was at 8 cm from the inter¬ 
action point. The green light emitted by the Lanex was 
selected using a 546 nm band pass filter and recorded by 
a 12 bits CCD. When in use, the Lanex screen excluded 
the electron spectrometer. In addition to the electron 
diagnostics, a Thomson parabola was used to detect pro¬ 
tons emitted from the target rear along the normal to 
the surface, as in a previous experiment [19j. The cut-off 
energy of protons was used as a reference to optimize the 
focal position of the target. 

The grating targets were produced at Czech Technical 
University, Prague by heat embossing of Mylar™ foils 
using a metallic master. The target material was cho¬ 
sen considering its high damage threshold for prepulses. 
In the following we show results obtained with targets 
having a resonant angle of ~ 30°, i.e. = 2A having 
assumed uo p uj. The average thickness was 10 pm and 
the peak-to-valley depth of the grooves 0.25 pm. Flat 
foils with the same average thickness were used for com¬ 
parison. In a limited number of shots also gratings with a 
resonance angle of 15° (X g = 1.35A) and 45° (X g = 3.41 A) 
were used, obtaining similar results. 

The electron emission from the front face of the tar¬ 
get changes dramatically between gratings irradiated at 
angles near resonance and flat foils. Fig. [2] shows the 
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FIG. 2: Images on the Lanex screen for simple flat target (top) 
and for grating target (bottom), both irradiated at pi — 30° 
incidence. A 3mm A1 foil was placed in front of the screen to 
filter out electrons with energy E < 1.7 MeV. The parabolic 
dashed lines give the local 6 angle corresponding to the posi¬ 
tion on the screen. 


spatial distribution of the electrons, as collected by the 
Lanex screen. The emission from the flat foil is rather 
diffused, with a “hole” in correspondence of the specular 
reflection of the pulse, as if electrons were swept away 
by intense light. The signal is maximum in an annular 
region around the hole. In contrast, for a grating at res¬ 
onance the emission is strongly localised on the plane of 
incidence (6 ~ 0°). The maximum intensity is close to 
the target tangent and is ~ 10 times larger than what 
observed for flat targets at the same angle of incidence. 
Two minima (holes in the image) are observed in the di¬ 
rections of specular reflection and first-order diffraction 
of the laser pulse (evidence of grating diffraction of the 
high-intensity pulse, which confirms the survival of the 
grating during the interaction, was also found in previ¬ 
ous measurements (HI). Local bending of the target or 
non-exact perpendicularity of the grating grooves to the 
plane of incidence may result in shot-to-shot fluctuations 
of the direction of maximum emission. Depending on 
the individual foil, the average angular shift in 6 was in 
the 1° — 5° range. An optimization of target design and 
alignment is foreseen to eliminate the fluctuations. 

The energy spectra were obtained placing the spec¬ 
trometer at p spec = 2°. The angle of incidence of the 
laser (pi was varied from 20° to 52°. Fig. [3] shows spectra 
obtained for (pi > 30° as for smaller angles no clear signal 
above the noise level was collected. The aforementioned 
fluctuations of the direction of the electron beam lead to 
a shot-to-shot variability of the intensity of the signal. 
Nevertheless, the most intense signals are detected only 
close to the resonance angle (^ 30°). Moreover, spectra 
collected at 30° and 35° are characterized by higher max¬ 
imum energies and a peculiar distribution with a dip at 
lower energies (3-4 MeV) and a broad peak at 5-8 MeV. 
On the higher energy side electrons with energy up to 
~ 20 MeV are detected. Whichever the angle of inci¬ 
dence, we never observed an electron spectrum above the 
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FIG. 3: Electron spectra collected with the detector at 2° 
from tangent direction, for several pulse incidence angles 
(from 30° to 52°). In the upper panels, all the collected shots 
are shown and a few of them are highlighted in color, while 
the others are plotted in light grey to show the shot to shot 
variations. 


noise level from flat targets at p spec ^2°. We also anal¬ 
ysed the electron spectra obtained irradiating the grating 
at (pi = 30° and changing the position of the spectrome¬ 
ter in the p spec = 1° — 35° range. Despite the shot-to-shot 
fluctuation, the shape of the spectra remained similar in 
all positions of the spectrometer (p spec < 20°,with a peak 
at 5-8 MeV. On the other hand, the intensity of the sig¬ 
nal monotonically decreased with respect to <p spec and the 
signal was visible up to <p spec — 30°, in agreement with 
the signal collected on the Lanex screen. 

Massively parallel 3D simulations were performed con¬ 
sidering flat targets irradiated at (pi = 30° and gratings 
irradiated at (pi = 30°, 35°, 40°. For reasons of compu¬ 
tational feasibility, the thickness £ t of the simulated tar¬ 
get was £ t = 1A and the electron density was n e = 50n c 
(where n c = 7rm e c 2 /e 2 A 2 is the cut-off density), while the 
other parameters corresponded to the experimental ones. 
The simulations were performed on 16384 cores of the 
FERMI supercomputer using the open-source particle- 
in-cell (PIC) code “PICCANTE” [A EH. The numeri- 
cal box size was 80A x 80A x 60A, wide enough for the 
boundaries not to affect the results. A spatial resolution 
of 70, 51, and 34 points per A in each direction and 50 
particles per cell were used. 

Fig. 0 shows a snapshot of a 3D simulation of a grat¬ 
ing target irradiated at its resonance angle (the inter¬ 
action geometry is shown in FiglTj) . The magnetic field 
component B z is represented together with the isosurface 
corresponding to the electron density of the grating tar¬ 
get. A surface wave propagating along the — y direction 
is excited and electrons trapped in the SP are accelerated 
along the target. Fig. [5] a) shows the simulated electron 
energy spectra dN/dE at p spe c = 2° for the flat target 
irradiated at (pi = 30° and grating targets irradiated at 
(pi = 30°, 35°, 40°. The signal for the flat target case is 
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FIG. 4: 3D particle-in-cell simulation snapshot of laser¬ 
grating interaction, showing the magnetic field component B z 
and an isosurface highlighting the target at time t = 45\/c 
after the beginning of the interaction. Only the z < 0 region 
of simulation box is shown in order to highlight the held dis¬ 
tribution. In the upper part of the picture, the distribution 
of B z is related to the diffracted pulse, while that in the lower 
part corresponds to a surface plasmon propagating in the — y 
direction. 



0 5 10 15 20 25 0 5 10 15 20 25 



0 10 20 30 40 50 60 70 80 90 


FIG. 5: Analysis of 3D simulations at time t = 45A/c after the 
beginning of the interaction, a) Simulated spectrum dN/dE 
at 4> sp ec = 2° for hat (F) targets and for gratings (G) at 
different values of the incidence angle fa. Particles with x < 
—0.125 (the average surface position of the target), \<p — 2| > 
0.5° or |0| > 0.5° are hltered out. b) Comparison of spectra 
between 2D and 3D simulation for a grating target irradiated 
at fa a? 30°. c) Electron angular distribution d 2 N/d(j)d0 for 
grating irradiated at 30°. Particles with x < —0.125 or kinetic 
energy Ek < 1 MeV are hltered out. 


very weak compared to grating spectra and it exhibits an 
energy cut-off which is ^ 10 x lower. The energy spec¬ 
trum for the grating irradiated at fa = 30° shows the 
peculiar spectral shape observed for fa = 30°, 35° in the 
experimental results (see Fig. G2 1, while for higher an¬ 
gles of incidence the low energy dip is not observed as 


for (pi >40° in the experiment. Fig. 0 b) compares the 
spectra obtained in 2D and 3D simulations. Only the 
3D simulation fully reproduces details of the spectrum 
such as the broad peak with low energy dip. Fig. [5] c) 
shows the simulated angular distribution on the screen, 
which also reproduces fairly well the experimental data 
(Figj2j), including the hole in correspondence of the angle 
of specular reflection. 

The 3D simulation also shows a correlation between 
electron energy and emission angle. Electrons at energies 
lower than the peak value are emitted at some angle with 
respect to SP propagation direction, so that integrating 
over the whole range of 6 the spectrum resembles the 
one observed in the 2D case. This is consistent with 
interpreting the fluctuations in the energy spectra (Fig.3) 
as related to those in the electron beam direction. 

In the simulation, a SP ~ 1 showing that the SP is rel¬ 
ativistic. Inserting such value and a = 50 in Eq.(j2]) we 
obtain £f — 25 MeV, <p e ~ 8° and cx 8A, in fair agree¬ 
ment with the observations. The linear scaling of £f with 
ao is also apparent in parametric 2D simulations. Details 
of the spectrum will be dependent on the distribution of 
injection velocities. In the supplemental material, movies 
from 2D simulations showing the acceleration of electrons 
moving along the surface are reported. In a very recent 
work which explores electron accerelation regimes in sur¬ 
face plasma wave [Hj, possible self-injection and phase¬ 
locking of electrons at relativistic intensity in the surface 
plasma wave is shown with a test particle approach. 

Our results may have application as an intense ultra- 
short electron source in the multi-MeV range, with char¬ 
acteristics not easily attainable with other techniques. 
Electrons with these parameters are potentially of in¬ 
terest for photoneutron generation [26H28] (aiming at 
extremely high fluxes) or ultra-fast electron diffraction 
( 29 I-EH (for imaging of ultrafast processes with electron 
diffraction). Besides these possible applications our re¬ 
sults confirm the possibility to extend the study and ap¬ 
plications of plasmonics into the high field, relativistic 
regime. 
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